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Abstract 

We present the Monte Carlo program BHAGEN95, for calculating the cross- 
section of the Bhabha scattering process at LEPl/SLC and LEP2 energies, usable 
with continuity from smah to large- angle configurations. We discuss some improve- 
ments in the event generator BHAGEN94, which is now part of the new code. 
In particular the weak and QCD corrections are implemented up to two loops for 
the relevant contributions, and the emission of one hard photon is treated exactly. 
We have included all the radiative corrections which are necessary to obtain, for 
a typical experimental event selection, a precision of 0.1-0.2% at small-angle. At 
large-angle we estimate a precision of 0.5%, with the exception of the region where 
the beam energy is a few GeV above the Z boson resonance, where it is up to 1%. A 
detailed comparison with other codes for both small-angle and large-angle Bhabha 
scattering is performed. 



PACS 12.15.-y Electroweak interactions 
PACS 12.15.Lk Electroweak radiative corrections 
PACS 12.20.-m Quantum electrodynamics 
PACS 12.20.Ds Specific calculations 
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1 Introduction 



Bhabha scattering is measured with remarkable precision at LEPl/SLC energies in two 
kinematical regions identified by the scattering angle: at small-angle for the luminosity 
monitoring and at large-angle for the measurements of the Z properties. 

Such results require a continuous effort on the theoretical side to maintain the the- 
oretical error on the prediction within a comparable precision and to provide a tool for 
realistic comparisons with the experimental data (namely an event generator). 

We have already given in [|T| the general description of our approach to the problem, 
while the features of the first version of our event generator BHAGEN are described in 
0, with instructions on how to get the program and how to use it. In we have con- 
sidered the higher order QED corrections, including the most relevant contributions and 
estimating the omitted terms; in a new version of our event generator, BHAGEN94, 
is also announced. 

Several improvements have been made subsequently to BHAGEN94, in order to in- 
crease and control its precision, and in some results for small-angle Bhabha scattering 
cross-section are presented with the label BG94-NEW and compared with the results from 
few other available computer programs: OLDBIS §, LUMLOG §, BHLUMI 2.01 [0] 
and 40THIEVES [Q. The conclusion from the comparisons with the above programs, 
and also from the use of the new event generator BHAGEN-IPH, which contains the com- 
plete hard photon matrix element 0, is that the hard photon part (even in the improved 
version of BHAGEN94) deserves a more accurate treatment to achieve the 0.1% precision, 
presently aimed at small-angle for luminosity monitoring. To obtain such a result a new 
procedure and the code BHAGEN95 are developed. BHAGEN95 (available on request 
from the authors) is actually a collection of three programs and in this paper we discuss 
in detail two of them, as the third one is already discussed separately [^. The current 
precision for the cross-section is heavily dependent on the chosen event selection and for 
our code can be estimated to be about 0.1%-0.2% for typical experimental event selection 
at small-angle for LEPl/SLC and LEP2 energies. In references [|lOl and |ll| our results 
are presented in a variety of different event selections and conditions (some of those pre- 
liminary results are updated in this discussion), and compared with the other presently 
available programs of similar precision: BHLUMI 4.03 0, OLDBIS+LUMLOG 1T|], 
SABSPV NLLBHA ig. 



The precision at large-angle depends on the beam energy and also on the event selection 
used; in the worst case (beam energy a few GeV above the Z peak) can be safely estimated 
to be about 1%, mainly due to the raise of the incertitude coming from the second order 
leading logarithmic correction in the s — t-interference terms, which are still unknown and 
usually are assumed to be the same as in the s-channel. For the other energy ranges at 
LEPl/SLC, and at LEP2, one can expect with relatively loose cuts a typical accuracy 
of 0.5%, which can be usefully compared with the other presently available programs, 
in the conditions they are suited for: ALIBABA 0, ZFITTER [|1|], TOPAZO [19 
BHAGENE3 |2g], UNIBAB [H, BHWIDE |2[, SABSPV pf. 



This paper goes through the following sections: implementation of an improved en- 
ergy distribution for the final photon emission in the resummed program BHAGEN94, 
in section 2; construction of the BG94-FO code, the 0{a) version of the event genera- 
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tor BHAGEN94, in section 3; implementation of necessary weak and QCD corrections 
in BHAGEN94 up to two loops, in section 4; construction of the BHAGEN95 code, in 
section 5; comparisons at small-angle, in section 6; comparisons at large-angle, in section 
7; conclusions, in section 8. 



2 Improved energy distribution for the final photon 
emission in the resummed program BHAGEN94. 

In the same notation as in [H, and p3|, the integrated cross-section is 



aie e 



e+e") 



AC- 



dQ. 



1) 



where the dressed differential cross-section in our approach is 
daiEf„e.) 



dQ. 
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Ef, is the beam energy, 9± the scattering-angle of the electron (positron), s = 4£'| and 
t± = — ssin^^. In the laboratory system the electron (positron) energy after initial 
radiation is E'^, after the scattering is E±, after the final radiation is E'^, and Emm is the 
minimum energy for the leptons. 

In terms of the variables in the LM-c.m.s. (the c.m.s. after initial radiation emitted 
collinearly to the initial direction, as introduced in [Q]), the incoming fermion energy 
E* = ^E^El, the scattering angle 6'*(E° , E^, 6^) and the jacobian J^) , the differential 
cross-section from Feynman diagrams, 0{a) complete and up to 0{a^) leading logarithm, 
is written as 



da{El,E^_,e_ 

dnl 



'dQ* 

dnl 



^daf{E\e* 

i=l 



dQ* 
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i + 4^(E5,r)l \i + c'^\E*,e*)] \i + s^^\E*,e* 



where iff^ ^ are the Born differential cross-sections for the various channels with 

vacuum polarization insertions (Dyson resummed as in Eq.(2.6) of i), 5^S{E,e) (defined 
in Eq.(16) of 0|) are the leading logarithmic corrections up to (9(a^), and C^^\Eb,9-) 
are the 0{a) complete corrections after subtraction of the terms already included in the 
previous corrections. 
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The function Q{cuts) accounts for the rejection procedure, according to the requested 
cuts on energies and angles. 



As usual, 



(3{E,e) 



j3,{E)+l3Ud) 
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with P{Eh,9±) ~ (3e{y/s/2) at large-angle and /3{Eb,6±) ~ l3eiy/~t±/2) at small-angle. 

The emission functions for initial (F") and final {F') emission in the notation of | 
and using the results of |2^, are 

F^{Eb, E^, 9±) 
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The radiator functions for initial (-D°) and final {D') emission are 



E. 



\ i/3(Sb,e±)-i 



D\Eb,Ele^) = -P{EM[l-^ 



and 



E. 



, X i/3(£;^,e±)-i 



(9) 



note that the leading logarithmic corrections, up to 0[oP'\ have been included in the 
factors + of the expression in Eq. (Bf). 

As the differential cross-section in Eq. does not depend on final lepton energies 
E'j^^ in our previous formulations the integrations on these variables were performed ana- 
lytically, assuming the further simplification that in the expression for D' in Eq. @ the 
values for P{E'^,9±) were taken constant at energy E*, i.e. P{E*,9±). 
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The simplification allows the analytical integration 

A^^^^*'''*^G(E*,r,A_) , (10) 

^i^^'''''''^G{E\e\^+) , (11) 

E 

where A± = 1 — Emin/ E± and the function G{E* ,9* ^ /S.±) is explicitly given in Eq. (23) 

of i. 

The procedure is rather good for small values of A±, but the extension to larger and 
more realistic values, corresponding to harder photon emission, causes a loss in precision, 
which justifies the difference with BHLUMI 2.01 and 40THIEVES [§ amounting up 
to 0.8%, as reported in for results under label BG94-OLD. 

In the new version of BHAGEN94 the (3 dependence on E'j^ is kept and the integration 
(or generation) is done numerically without further approximations. The new values 
(reported in |^ under label BG94-NEW) reduce the difference with the above mentioned 
programs to about 0.25%. 

The new approach imply the generation of two more variables for the final fermion 
energies E'j^, for the proper description of the peaks in Eq. (^ the variables y'^ are intro- 
duced 

y'. ^ (i - 1) . (12) 

To generate within a unit volume, the variables are changed to G (0, 1) 

y± = ^±[^-^) ■ (13) 

Note that, even if in the generation procedure the function (3 is still taken constant in 
energy, in the actual calculation the exact value is computed and used to determine the 
cross-section. 




3 Construction of the BG94-FO generator. 

We have extracted from BHAGEN94 a first-order event generator — i.e. including only 
corrections up to 0{a.) — that we call BG94-FO, keeping as much as possible unchanged 
the structure of the program. 

That is achieved by analytically expanding in a all the formulae and then keeping the 
0{q) terms only; the separation between soft and hard radiation is reintroduced through 
the small emitted energy fraction cut e, whose size has to be of the order of the accepted 
error, typically e < 10^^. To allow the expansion in a, the integrations in Eq. (^ are 
split into 

Eb Eb(l-e) Eb 

j dEl= I dEl + I dEi, (14) 



4 



and 

E± E±{l-e) E± 

J dE'^= J dE'^ + J dE'^. (15) 

All the resulting terms, which involve more than one integration in the hard photon 
regions (-Emm, Eb{l — e)) and {Emm, E±{1 — e)), are omitted, being of higher order. 

In the soft photon regions {Eb{l — e), Ef,) and {E±{l — e), E±) and for smooth functions, 
E\ can be approximated by Ei, and E'^ by E±, while special care has to be devoted to 
the peaking behaviour of the functions and D' . Moreover from the definition of E^ 
(Eq. (3.6) of 0) one has that for E\ ^ also E± ~ E^, and 6+ = 6^, so Eq. (|) 
becomes 
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where the first term accounts for soft emission only, the second and third account for hard 
emission from final electron or positron, the fourth and fifth for hard emission from initial 



5 



electron or positron. The functions and F' are now intended at 0{a), that is omitting 
the P term in curly braces of Eq. @ and Eq. (|). 

The analytical integration and subsequent expansion in a can be easily done for the 
initial radiation 



iPi'^'^M = l + \(3{EM\Tie + 0{a') , (17) 



and in the case of final radiation the result is similar, apart from corrections of order e 

^D'(E.. E'. . f)^) = fh(^(E'"(>i)^n(f] = 1 + ^ 



Eb{l-e) 



E 



^D'{E,,E'^,e±) = e5/5(^'"^±) + 0(e) = l + \l3{E,,e±)hie + 0{a') + 0{e) . (18) 



In Eq. d) the expansion in a is performed and only 0{a) terms are retained: the vacuum 
polarization contribution 6yp is defined in Eq.(2.3) of 6^^{Eb) is taken up to 0{a) 
from Eq.(16) of 0, and C^'''\Eb,6^), as indicated in 0, is the nonleading part of the 
0{a) correction C^'\Ef„e_), given explicitly in ||2^. When an 0{a) correction from 



above integrals occurs only the Born cross-section is retained, and if no longer necessary 
the transformation to LM-c.m.s. is released. 

With all that the cross-section at 0{a) becomes 
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with accordingly E* = J E^El and 6* = 9*{Eb, E^, 9-). 



In Eq. (^) the last two lines are due to hard initial emission, the second and third 
lines are due to hard final emission; in the approximation of constant (3 and symmetrical 
cuts, the hard final emission contributions can be integrated analytically giving 
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with A = 1 — Emin/ Eb. In Table 1 of 0] Eq. ( pi]) is used to produce the results presented 
under the label BG94-FO-OLD, which differ from the exact ones up to 0.87%, while the 
direct generation from Eq. (|TU|) gives the results under the label BG94-FO-NEW, which 
have only up to 0.28% difference from the exact ones; these are under label OLDBIS [|] 
and BG94-FO-EXACT (obtained using BHAGEN-IPH |] for the hard photon emission 
part) and are well in agreement inside the OLDBIS technical precision of 0.02%. 



4 Implementation of weak and QCD corrections. 

The previous versions of the program |[^, 0, including the latest BIIAGEN94 0, did 
not contain the Zee weak vertex correction, while the self energy corrections were included 
in an approximate form. For an accuracy better than 1% at large- angle those complete 

m 



corrections and also some higher order corrections |^ are necessary. After their 

inclusion the error due to the neglected weak one loop corrections is about 0.06% around 
the Z boson resonance (LEPl) and 0.1% at LEP2 energies, so smaller than the other left 
over photonic corrections. The use of the program at much higher energies (say 1 TeV) 
requires an update of the weak library for having a precision better than 2%. 

For the one loop Zee vertex corrections and self-energy corrections we use the results 
of ||25|, but we introduce the Dyson resummation of the self energy, as described for 
s-channel in |^ (we do actually the Dyson resummation in all channels). The higher 



orders corrections are included as outlined in [Q, using the available results for weak 
2§ and QCD corrections [gUf. 

For completeness we report here the formulae of these corrections as they come using 
our notation as in 111, 0|, P] and 
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while functions 5-'|,y(±), 5|f,|/(±) (i = 1,2,3), 11'^^ (t) and constants are defined in 
the second reference of in the same notation. As usual we take 



sin^(6'vF) = = 1 - c: 
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(31) 



8 



The above formulae take into account the Zee weak vertex corrections and mixing 
self energy contribution. To include completely the Z boson self energy corrections it is 
sufficient to change in our approach the two functions x(s) and The new ones are 

IQsfycfy 1 + (s) s - M| + isTz/Mz 



where Il^{s) = Re[S^(s)]/(s — M|) and is again defined in All the higher 

order corrections are included into the and IP^ functions, and through them into the 
formulae ( p3H30| ) and (p^-|33|), which remain unchanged. 

The value of is calculated using Eq. (|3ll) with Mz as an input parameter. In 
the code the value of M^y can be supplied by the user, or is calculated using the Fermi 
constant Gp as an independent parameter. In the latter case we implement the approach 



outlined in ||2^ , which allows to get the value of the W boson mass with higher accuracy 



than from its direct measurement. For this better precision it is necessary to include 



the already mentioned weak corrections up to two loops to the Ap parameter ||2^ and 
also the QCD corrections to the self energy functions of the bosons (which are needed to 
properly calculate Ar ) |2^. Checking with the values of the W boson mass presented 
in |]2B[ as a function of the mass of the top quark (at fixed values of rriH = 300 GeV 
and as{Mz) = 0.125), we obtain a very good agreement around mt = 175 GeV (only a 
small difference (0.03%) is appreciable above the unrealistic value mt = 225 GeV). 

Actually the inclusion in BHAGEN94 (and as a consequence also in BHAGEN95) 
of the higher order QCD corrections is done in an approximate form. In particular for 
the first two fermion generations (u,d,s and c quarks, for which the limit s, —t ^ is 
valid, with m the quark mass) we use the multiplicative factor (1 + ^ + ...) in front of 
the contribution due to one quark loop in the self energy functions of the bosons. For 



the t,b doublet we use the approximated formulae, presented in p9[, for the doublet 
contribution to Ar in the Mw calculation. On the contrary in the Bhabha cross-section 
calculations is necessary the knowledge of the separate contributions 11^ and 11'^^ to the 
boson self energies, we use them in the form described in |^6| (BHM/WOH approach), 
i.e. we include the leading Ap^'^ corrections. In this way the corrections are calculated 
in the point s = M| and their s, t dependance (which is known however only in the 
QCD two loop contribution) is not taken into account. Also some known nonleading 
corrections are not included. All this considered the approximation is sufficient to reach 
at large-angle the projected accuracy of 0.5% in the cross-section calculations. Moreover 
other not included QED two loop corrections are expected to be much bigger than the 
left over QCD corrections, like the next to leading correction, which is not yet calculated, 
or the inferred leading term in s — t channel interference, assumed to be identical to the 
s-channel. 
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5 The BHAGEN95 code. 

BHAGEN95 is a collection of three programs to calculate the cross-section for Bhabha 
scattering from small to large scattering angles at LEPl and LEP2 energies. In its present 
form the integrated cross-section cr(BHAGEN95) for a given selection of cuts is calculated 

as 

cr(BHAGEN95) = ct(BHAGEN94) - (t^(BH94-FO) + (t^(BHAGEN-IPH) . (34) 

cr (BHAGEN94) is the integrated cross-section, based on leading exponentiated 
formulae, obtained with the Monte Carlo event generator BHAGEN94, with the imple- 
mentations described in the previous sections 2 and 4. The use of collinear kinematics of 
initial and final radiation leaves some approximation in the angular distribution, which 
limits the accuracy particularly in the region of hard photon emission, as remarked at the 
end of section 2. 

{BH.94-FO) is the integrated cross-section of 0{a) for one hard photon emission, 
obtained selecting the appropriate part in the Monte Carlo event generator BH94-FO, the 
0{a) expansion of BHAGEN94 described in section 3. 

(T^(BHAGEN-IPH) is the integrated cross-section obtained with the one hard photon 
complete matrix element and exact kinematics, implemented in the Monte Carlo event 
generator BHAGEN-IPH 0. 

The subtraction of a^(BH94-FO) and its substitution with (t^(BHAGEN-IPH) is 
performed to reduce the error in the cross-section, coming from the approximation in the 
contribution due to the one hard photon emission in BHAGEN94, as discussed in . In 
this way the one hard photon emission is treated exactly in BIIAGEN95. 

The choice of the energy threshold e (in units of the beam energy) for a photon 
to be considered hard should not affect the difference between (BHAGEN-IPH) and 
(T^(BH94-FO). However the approximate hard photon treatment in BH94-FO implies 
that a correction of the order of e is missing, so it is necessary to choose e well below 
the desired accuracy, typically in the range of 10~^ — 10~^. But even so the difference 
is really independent on e only if the soft photon treatment in BH94-FO is the same as 
in the exact program BHAGEN-IPH. This is indeed the case. In fact changing e from 
10~^ to 10~^ the separate cross-sections are growing about 20-30% for the event selections 
described in [^, while the difference is stable within statistical errors. 



The three programs provide cross-sections, which are summed as in Eq. (|^) or used to 
obtain other quantities, such as forward-backward asymmetry. Although the constituent 
programs are separately genuine event generators, the discussed combination can be used 
for Monte Carlo integration only. 

At small-angle we estimate the accuracy in the cross-section evaluation, due to the 
uncontrolled higher orders terms 0{a'^L) and 0{a^L^) and to the incertitude in (9(a^L^) 
s — t interference, to amount comprehensively to about 0.1%. The further error, due to 
the approximate two hard photon contribution (strongly dependent on the imposed cuts) 
is estimated on the basis of the calculated correction for the one hard photon contribution 
times (3e{s) — 0.1, to account for the increase in perturbative order. 
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In Table |^ are presented the values of the ratio R, defined as 

^^^ a^(BHAGEN-lPH) - o-^(BH94-FO) 
^ = a(BHAGEN95) ' ^^'^ 

which is the correction (in %) introduced in BHAGEN95, to account for the approximation 
in the one hard photon contribution in BHAGEN94. The correction is given for the event 
selections used in ||T0[ , in the same notation as for the t-channel comparison. As already 
said, for /3e(s) ~ 0.1 the same numbers in Table |l| are an estimation (in per mil) of the 
inaccuracy, due to the approximate treatment of two hard photon emission. The values 
in Table |l] are acceptably small for cuts closer to the ones in experiments, which are 
of the angular asymmetric WN type, with 0.3 < Zmin < 0.7, where Zmin is the energy 
threshold for the final clusters for accepting or rejecting the events. Indeed in such cases 
the (included) one hard photon corrections are found to be below 1.4 % at LEPl, so that 
the corrections to two hard photon emissions are expected to be below 1.4 per mil. For 
the values at Zmm = 0.9 and calorimetric event selection, these corrections are expected 
to be much bigger (at a few per mil level). 

All included we estimate at small-angle an accuracy of the order of 0.1%-0.2% for 
typical experimental cuts for both LEPl and LEP2 energies. 

At large-angle we estimate the accuracy of the (9(a^L^) s — t interference contribution 
up to 1% (depending on energy and cuts) at LEPl, but much smaller at LEP2. The 
error, coming from the approximate treatment of two hard photon emission, is estimated 
as explained above for the small-angle case, and is smaller for more stringent acoUinearity 
cut. All included we estimate an accuracy of the order of 1% in the worst case at LEPl, 
when the beam energy is a few GeV above the Z boson peak, while typical accuracy for 
relatively loose cuts is 0.5%. 

The three programs run separately. They provide initialization and fiducial volume 
definition according to input parameters, then starts the generation of events according 
to appropriate variables, which smooth the cross-section behavior. Rejection is performed 
through the routine TRIGGER, where the special cuts can be implemented. The programs 
stop when the requested number of accepted events is reached or alternatively when the 
requested accuracy is obtained. 

The following data have to be provided in input: mass of the Z, mass of the top quark, 
mass of the Higgs, value of as{Mz), value of Tz, the beam energy Ebeam, the minimum 
energy for final leptons Emin (larger than 1 GeV), minimum and maximum angle for 
the scattered electron (positron) with the initial electron (positron) direction, maximum 
acoUinearity allowed between final electron and positron, number of accepted events to be 
produced, numbers to initialize the random number generator. The following possibilities 
are also provided: i) to switch on or off the leading contribution from virtual and soft 
emitted pairs , ii) to calculate separately the different channel contributions (useful for 
tests), iii) the recording of the generated events of each component program in a separate 
file. 

For 0{a) programs the minimum and maximum energy allowed for the photon has to 
be specified. The input of BHAGEN-IPH requires also the maximum acoplanarity, and 
minimum angles of the emitted photon with initial and final fermion directions, if the 
contributions with the coUinear photons are to be excluded. 
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Each program returns the input parameters and the values of the cross-section obtained 
with weighted and unweighted events, with the relative statistical variance (one standard 
deviation). Of course, due to the efficiency, the weighted cross-section is usually much 
more precise than the unweighted one. The total integrated cross-section is then calculated 
according to Eq. (13^ 





0.1 


0.3 


0.5 


0.7 


0.9 


BAREl(WW) 




-0.24 


-0.29 


-0.23 


-0.11 


CALOl(WW) 




-0.46 


-0.75 


-1.15 


-2.29 


CAL02(WW) 


-0.38 


-0.60 


-1.01 


-1.66 


-3.57 


SICAL2(WW) 


+0.35 


+0.26 


+0.005 


-1.02 


-3.34 


CAL02(NN) 


-0.54 


-0.73 


-1.06 


-1.76 


-3.66 


CAL02(NW) 


+0.01 


-0.22 


-0.64 


-1.39 


-3.40 


CAL02(WW-all inch) 


-0.36 


-0.57 


-0.96 


-1.59 


-3.42 


CAL02(NN-all inch) 


-0.52 


-0.70 


-1.02 


-1.69 


-3.51 


CAL02(WN-all inch) 


+0.00 


-0.21 


-0.61 


-1.33 


-3.26 


CAL03(NN) 


-1.14 


-1.29 


-1.56 


-2.19 


-4.38 


CAL03(NW) 


-0.21 


-0.41 


-0.77 


-1.61 


-3.95 


CAL02(WW) 




-0.61 


-1.02 


-1.68 


-3.62 


SICAL2(WW) 




+0.25 


+0.04 


-1.03 


-3.40 



Table 1: Correction R, defined in Eq. (|35D, due to the approximate treatment of the one hard 
photon emission in BHAGEN94, in % of the BHAGEN95 cross-section. All event selections are 



defined in |10| in the same notation. First 9 rows are for LEPl and last 4 rows are for LEP2 
energies. The option 'all inch' means that all possible contributions are included, while the 
remaining tests are done only for t-channel contributions with vacuum polarization switched off. 



6 Comparisons at small-angle. 

The 0{a) results are tested with other programs in for bare event selection and in 
|T0[] also for calorimetric event selection, where the precision of 0.03 % is confirmed, as a 
consequence of the agreement of different programs within statistical errors. 

In 1^ the results of BHAGEN94, including the improvement outlined here in section 
2, are presented for a comparison of the cross-section with the final electron and positron 
scattering angles d± in the range 3° < ^± < 8°. From this comparison it is clear that we 
cannot reach a precision better than 0.3% using only the structure function approach. 

To have a better precision the code BHAGEN95 is settled, with the features described 
in the previous section. Several preliminary results are already presented and compared 
in JlOl , so we correct the few which are revised, but we do not repeat here the unchanged 
ones. 

In Fig. 16 of ||10|, in the comparison at small-angle (t-channel only) the bare event 



selection case (BAREl, WW) shows good agreement (almost inside the 0.1% box) up to 
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Zmin = 0.9 included among all the programs (BHAGEN95, BHLUMI 4.03 |g, OLD 



BIS+LUMLOG [|T5l, SABSPV 0, NLLBHA |Tj). 

For calorimetric event selections BHAGEN95 and OLDBIS+LUMLOG have very close 
results for every event selection cut, but they agree (within 0.1% at LEPl and 0.2% at 
LEP2) with BHLUMI and SABSPV only for values of the parameters close to the real 
experimental ones (NW or WW, 0.3 < Zmin < 0.7). The difference is larger for more 
severe cuts (as in the NN case) and for Zmin = 0.9. 

It is proposed in [|1^] that for OLDBIS+LUMLOG the difference is due to the so 
called 'classical limit' (i.e. zero radiation emission limit) which could be different in 
higher orders. 

For BHAGEN95 we believe that the difference is due to the approximate treatment of 
two hard photons illustrated in previous section, while the 'classical limit' is sufficiently 
accurate. In fact for calorimetric event selection, even for Zmin 1-0, that limit is not 
approached, as almost collinear hard photons can join the final lepton in the cluster. The 
'classical limit' is approached only for BAREl event selection, where all the programs are 
in substantial agreement even at Zmin = 0.9, and remarkably there the estimated error in 
BHAGEN95 on the basis of Table |I| is only about 0.01 %. On the contrary for calorimetric 
event selection and Zmin = 0.9 hard photons are included within very stringent cuts on 
phase space, as the cluster opening is very narrow. The structure function approach (used 
in the part of BHAGEN95 called BHAGEN94) for this configuration is expected to be 
inaccurate, as it is not able to mimic very sophisticated cuts. 

In conclusion, at small-angle for the test cross-section called 't-channel only', the 
results of BHAGEN95 are in very good agreement with those of OLDBIS+LUMLOG, and 
inside the errors also with those of BHLUMI and SABSPV. The accuracy of BHAGEN95 
is very much dependent on the event selection used, and is of the order of 0.1% - 0.2% for 
typical experimental cuts, while for more stringent cuts it can amount up to say 0.5%, 
with the source of the error well understood. 

The results for the complete cross-section of BHAGEN95, with all the other contribu- 
tions included, are presented and compared in Table 18 and in Fig. 17 of [|l^, repeated 
here in Fig. for completeness Q. 

It is interesting to compare the difference Aa between the complete results (all in- 
cluded) of Table 18 and the 't-channel only' results of Tables 14 and 16 of |jlO| for 
BHLUMI and BHAGEN95 (for SABSPV values the comparison has no significance due 
to larger statistical errors). In Table |^ are presented the values for different angular ranges 
(WW, NN, WN) and energy threshold taken from Table 18 of ^ for BHLUMI 
and calculated as indicated above for BHAGEN95. 

The difference Act is mainly due to vacuum polarization correction, presumably im- 
plemented in the same way in both programs, using also the parameterization in |30 



and to the Z — 7 interference term. This latter contribution at small-angle is in the range 
of 1 per mil of the total cross-section, depending mostly on the angular range allowed and 
much less on the other details of the selection. 



^The values of BHAGEN95 in Fig. 2 of ||T^ differ slightly from the ones in due to a second bug 
introduced in the rush-fixing of a previous (actually inactive) bug in the routine TRIGGER. The corrected 
results confirm those already presented in [10] . 
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In the last column of Table |^ is presented in per mil the variation 

A(t(BHAGEN95) - Aa(BHLUMI) 



V = 1000- 



a(BHLUMI,all incl^ 



(36) 



between the complete cross-section obtained with BHLUMI and the one obtained with 
BHAGEN95, due to the difference in vacuum polarization implementation and Z — 7 
interference correction. It can be seen that in the experimentally interesting region the 
variation is statistically significant, always positive and increasing with Zmin, and it is 
almost up to 0.1%. 





A(t(BHLUMI,WW) 


Aa(BHAGEN95,WW) 


V(WW) 


0.1 


5.140(8) 


5.203(11) 


0.46(10) 


0.3 


5.126(8) 


5.197(14) 


0.52(12) 


0.5 


5.100(8) 


5.195(14) 


0.70(12) 


0.7 


4.994(8) 


5.125(16) 


0.99(14) 


0.9 


4.627(8) 


4.821(16) 


1.57(15) 




A(t(BHLUMI,NN) 


Aa(BHAGEN95,NN) 


V(NN) 


0.1 


3.751(7) 


3.802(14) 


0.51(16) 


0.3 


3.742(7) 


3.801(14) 


0.60(16) 


0.5 


3.728(7) 


3.799(14) 


0.72(16) 


0.7 


3.678(7) 


3.774(14) 


0.99(16) 


0.9 


3.430(7) 


3.579(14) 


1.64(18) 




Aa(BHLUMI,WN) 


A(t(BHAGEN95,WN) 


V(WN) 


0.1 


3.883(4) 


3.920(13) 


0.36(14) 


0.3 


3.873(4) 


3.919(13) 


0.45(14) 


0.5 


3.854(4) 


3.916(13) 


0.61(14) 


0.7 


3.779(4) 


3.869(13) 


0.90(14) 


0.9 


3.478(4) 


3.624(14) 


1.59(15) 



Table 2: The difference Ao", in nb, between the complete cross-section and the 't-channel only' 
contribution, for BHLUMI and BHAGEN95, for several event selections (WW, NN, and WN 
angular range), for 2Ei, = 92.3 GeV. In the last column is the variation V, defined in Eq. (|3^), 
giving in per mil the difference due to vacuum polarization implementation and Z— 7 interference 
correction of BHAGEN95 respect to BHLUMI. In brackets is the statistical error on the last 
digits. 



It seems difficult that the difference comes from vacuum polarization implementation, 
so we have investigated the Z — 7 interference contribution. 

In the next section are reported comparisons at large-angle: BHAGEN95 is in a good 
agreement (about 0.5%) with the other programs (notably with ALIBABA for BARE 
event selection), where the contribution of the Z — 7 interference term amounts up to half 
of the cross-section. 

The accuracy of the Z 



7 interference term included in BHLUMI is tested in 31 
relaying mostly on the comparison with the ALIBABA code, so we have obtained results 
in the same conditions for BIIAGEN95 for comparison. The results of ALIBABA and 
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the Born cross-sections, to which the resuhs are normahzed, were taken from for 
old generation luminometers (LCAL: larger angular range between 3.3° and 6.3°) and 
from for new generation luminometers (SCAL: smaller angular range between 1.5° 
and 3.15°, more similar to the ones used in |1^), while the results of BHLUMI are taken 
from 



in the same conditions. In Table ^ the results of ALIBABA, BHLUMI and 
BHAGEN95 in the SCAL angular range are inside the 0.01% difference, and in the LCAL 
angular range are inside the 0.05%. This test is done for very loose cuts: no cut on photon 
energy and angles is applied and only the angular ranges of leptons are restricted. In the 
results of Table ^ the agreement between BHLUMI and BHAGEN95 is much better for 
the looser cuts, than for the more severe ones, and this could be the reason for the better 
agreement in Table ^ . 

It is however difficult to identify with certainty the source of the difference without a 
more detailed analysis, particularly if it comes from a different use of the same electroweak 
parameters, but it would be necessary in case the accuracy has to be better than 0.1%. 





LCAL angular range 


SCAL angular range 




BHAGEN95 


ALIBABA 


BHLUMI 


BHAGEN95 


ALIBABA 


BHLUMI 


89.661 


0.757 


0.778 


0.794 


0.165 


0.172 


0.175 


90.036 


0.779 


0.799 


0.816 


0.169 


0.177 


0.180 


90.411 


0.724 


0.747 


0.754 


0.156 


0.164 


0.165 


90.786 


0.527 


0.545 


0.538 


0.112 


0.117 


0.116 


91.161 


0.175 


0.187 


0.158 


0.035 


0.036 


0.032 


91.563 


-0.206 


-0.206 


-0.247 


-0.048 


-0.050 


-0.057 


91.911 


-0.473 


-0.479 


-0.521 


-0.105 


-0.110 


-0.116 


92.286 


-0.602 


-0.609 


-0.647 


-0.132 


-0.136 


-0.143 


92.661 


-0.642 


-0.650 


-0.678 


-0.141 


-0.145 


-0.150 



Table 3: Z — -y interference term (in % of Born cross-section) for BHAGEN95, ALIBABA and 
BHLUMI for two different angular ranges (LCAL and SCAL) as a function of the beam energy. 



7 Comparisons at large-angle. 



At large-angle some comparisons were presented for LEPl and LEP2 energies in pO[ . 
However, as mentioned also there, the higher orders weak and QCD corrections were just 
implemented in BHAGEN95 and the results were very preliminary. Unfortunately some 
bugs were present in the program and, after their correction and test, the new results of 
BHAGEN95 can now be compared with the ones of the other programs ALIBABA |T7 

UNIBAB 



TOPAZO lH], BHAGENE3 |g, UNIBAB m, BHWIDE [||], SABSPV [0. 

We report here the new results of BHAGEN95 in Fig. |^ and Fig. ^ for LEPl and in 
Fig. ^ for LEP2 energies (corresponding respectively to Fig. 19, 20 and Fig. 21 in [|T^ 
in the same notation), for scattering angles 40° < ^_ < 140° and 0° < ^+ < 180°. 

We estimate that at large-angle the accuracy of the program BHAGEN95 is 0.5% 
everywhere, except in the region of a few GeV above the Z boson peak, where, due to 
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the absence of the calculation of the leading term in s — t channel interference, 

it can be assumed to be up to 1% f^. 

In Fig. 0, for LEPl energy and BARE type event selection the results of BHAGEN95 
are in agreement with most of the programs inside 0.5% for energies under and on top 
of the Z boson peak, while above the resonance and with more stringent acoUinearity 
cuts (10°) there is a noticeable agreement with ALIBABA, but some difference from 
other programs, which however is still less than 1%. The same pattern is kept for the 
calorimetric event selection CALO, shown in Fig. ^, where unfortunately results from 
ALIBABA are no longer available. 

In Fig. ^, for LEP2 energies and CALO event selection, there is a remarkable clustering 
in a 2% interval of most of the programs. 



8 Conclusions. 

The Monte Carlo integrator BHAGEN95 (available on request from the authors) calculates 
the cross-section of the Bhabha scattering with continuity from small to large-angle and 
from LEPl to LEP2 energies. 

The program contains all the necessary corrections to provide results, for typical ex- 
perimental event selections, with the precision required to usefully compare with experi- 
mental measurements and other theoretical results. We estimate a precision of 0.1-0.2% 
at small-angle, in the angular range of the new generation of LEP luminometers, mainly 
due to the inaccuracy in in the two hard photon emission contribution. At large-angle for 
both LEPl and LEP2 energies, we estimate a precision of 0.5% with the exception of the 
region of the beam energies a few GeV above the Z boson resonance, where the relevance 
of the still missing calculation of the leading (9(a^L^) radiative corrections to the Z — 7 
interference term spoils the accuracy up to 1%, as in every other present calculation. 

We have performed a detailed comparison of the photonic and weak corrections for 
both small-angle (angular range of the new generation of LEP luminometers) and large- 
angle Bhabha scattering. 

At small-angle we obtain an agreement better than 0.02% at 0{a) with the program 
OLDBIS for every type of event selection. In the comparison beyond the 0{a), but lim- 
ited to the t-channel only, with vacuum polarization and Z — 'j interference switched off, 
the results of BHAGEN95 agree with the ones of OLDBIS+LUMLOG better than 0.03%, 
for typical experimental event selection, and have a difference less than 0.1% with BH- 
LUMI and SABSPV in these conditions. For event selections less realistic for experiments 
the results of BHAGEN95 (and also of OLDBIS+LUMLOG) can differ from the ones of 
BHLUMI and SABSPV of a few per mil. We believe to understand the difference (at 
least in the BIIAGEN95 case) as due to the inaccuracy in the two hard photon emission 
contribution, which depends on the details of the event selection used and is smaller for 
looser cuts. When all possible contributions are switched on the agreement between BH- 
LUMI, SABSPV and BHAGEN95 is inside 0.1% for realistic experimental event selection, 
as presented in Fig. (same as in Fig. 17 of |l3l)- The better agreement in this case is 



due to a slightly larger (but less than 0.1%) contribution of BHAGEN95, coming mainly 
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from the implementation of the vacuum polarization and/or Z — 7 interference correc- 
tion. We have investigated this small difference, but it is impossible to disentangle the 
contributions without ad hoc runs of all the compared programs. 

At large-angle and in the LEPl energy range the corrected results of BHAGEN95 
are shown in Fig. |^ and in Fig. ^; the difference with the other results is within 0.5% 
around the Z boson resonance, being bigger for energies a few GeV above the peak, where 
BHAGEN95 is in agreement with ALIBABA within 0.5% and within 1% with other codes. 

For LEP2 energies although the differences between the codes are larger (up to 2% for 
the cluster of the results of TOPAZO, BHWIDE, SABSPV and BHAGEN95), as shown 
in Fig. ^, the precision is comparable with the foreseen experimental accuracy. 
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Figure 1: Monte Carlo results for various symmetric/asymmetric versions of the CAL02 ES, for 
matrix elements beyond first order. Z exchange, up-down interference and vacuum polarization are 
switched ON. The center of mass energy is y/s = 92.3 GeV. Not available x-sections are set to zero. 
In the plot, the 0{a'^)JJp^ cross section (Tg^^ ^''O'^ BHLUMI 4.x is used as a reference cross section. 
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No. 


EcM 


TOPAZO 


BHWIDE 


BHAGFNE3 


ALIBABA 


BHAGEN95 










( 


a) BARE acolmax — 10° 










1. 


88.45 


.4579 ± 


0003 


.4560 ± 


0004 


.4495 ± 


0016 


.4575 ± 


0003 


.4578 ± 


0002 


2. 


89.45 


.6452 ± 


0002 


.6429 ± 


0006 


.6334 ± 


0023 


.6440 ± 


0003 


.6445 ± 


0003 


3. 


90.20 


.9115 ± 


0002 


.9087 ± 


0008 


.8997 ± 


0033 


.9090 ± 


0004 


.9095 ± 


0004 


4. 


91.19 


1.1846 ± 


.0002 


1.1797 ± 


.0010 


1.1847 ± 


.0033 


1.1840 ± 


.0004 


1.1822 ± 


.0005 


5. 


91.30 


1.1639 ± 


.0002 


1.1592 ± 


.0009 


1.1667 ± 


.0033 


1.1636 ± 


.0005 


1.1619 ± 


.0005 


6. 


91.95 


.8738 ± 


0002 


.8711 ± 


0007 


.8856 ± 


0028 


.8769 ± 


0003 


.8742 ± 


0004 


i' . 


93.00 


.4771 ± 


0002 


.4761 ± 


0005 


.4808 ± 


0019 


.4814 ± 


0001 


.4796 ± 


0002 


8. 


93.70 


.3521 ± 


0002 


.3512 ± 


0004 


.3521 ± 


0013 


.3556 ± 


0001 


.3550 ± 


0001 




(h) BA] 


IE acol„r,x ~ 25" 


1 

1 . 


88.45 


.4854 ± 


0003 


.4808 ± 


0005 


.4699 ± 


0016 


.4833 ± 


0003 


.4826 ± 


0002 


2. 


89.45 


.6746 ± 


0003 


.6699 ± 


0006 


.6593 ± 


0023 


.6727 ± 


0003 


.6710 ± 


0003 


3. 


90.20 


.9438 ± 


0003 


.9387 ± 


0008 


.9279 ± 


0033 


.9425 ± 


0003 


.9384 ± 


0004 


4. 


91.19 


1.2198 ± 


.0003 


1.2130 ± 


.0010 


1.2169 ± 


.0034 


1.2187 ± 


.0004 


1.2133 ± 


.0005 


5. 


91.30 


1.1989 ± 


.0003 


1.1924 ± 


.0010 


1.1995 ± 


.0034 


1.1982 ± 


.0004 


1.1928 ± 


.0005 


6. 


91.95 


.9054 ± 


0002 


.9011 ± 


0007 


.9124 ± 


0026 


.9089 ± 


0003 


.9014 ± 


0003 


7. 


93.00 


.5040 ± 


0002 


.5013 ± 


0005 


.4996 ± 


0019 


.5054 ± 


0002 


.5027 ± 


0002 


8. 


93.70 


.3777 ± 


0002 


.3749 ± 


0004 


.3689 ± 


0013 


.3782 ± 


0001 


.3771 ± 


0001 
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Figure 2: Monte Carlo results for the BARE ES, for two values (10° and 25°) of acollinearity cut. 
Center of mass energies (in GeV) close to Z peak. In the plots, the cross section ctref from TOPAZO 
is used as a reference cross section. Cross sections in nb. 
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No. 


EcM 


TOPAZO 


BHWIDE 


RHAGENE3 


UNIBAB 


BHAGEN95 




(a) CALO acolmax — 10° 


1. 


88.45 


.4533 ± 


0004 


.4523 ± .0004 


.4467 ± 


0008 


.4490 ± 


0010 


.4524 ± 


0001 


2. 


89.45 


.6387 ± 


0004 


.6377 ± .0006 


.6302 ± 


0011 


.6358 ± 


0012 


.6370 ± 


0002 


3. 


90.20 


.9023 ± 


0003 


.9016 ± .0008 


.8920 ± 


0015 


.9021 ± 


0014 


.8990 ± 


0003 


4. 


91.19 


1.1725 ± 


.0001 


1.1707 ± .0010 


1.1767 ± 


.0021 


1.1772 ± 


.0016 


1.1689 ± 


.0004 


5. 


91.30 


1.1520 ± 


.0001 


1.1505 ± .0009 


1.1571 ± 


.0020 


1.1559 ± 


.0016 


1.1491 ± 


.0004 


6. 


91.95 


.8649 ± 


0001 


.8646 ± .0007 


.8795 ± 


0015 


.8689 ± 


0012 


.8660 ± 


0003 


i' . 


93.00 


.4723 ± 


0001 


.4725 ± .0005 


.4796 ± 


0008 


.4733 ± 


0008 


.4761 ± 


0001 


8. 


93.70 


.3486 ± 


0001 


.3486 ± .0004 


.3507 ± 


0006 


.3486 ± 


0007 


.3526 ± 


0001 




(b) CALO aco/„,„., ^ 25" 


1 

1 . 


88.45 


.4769 ± 


0004 


.4742 ± .0004 


.4696 ± 


0008 


.4733 ± 


0010 


.4741 ± 


0001 


2. 


89.45 


.6638 ± 


0003 


.6615 ± .0006 


.6556 ± 


0011 


.6619 ± 


0012 


.6599 ± 


0002 


3. 


90.20 


.9297 ± 


0003 


.9278 ± .0008 


.9207 ± 


0012 


.9302 ± 


0014 


.9236 ± 


0003 


4. 


91.19 


1.2025 ± 


.0003 


1.1994 ± .0010 


1.2074 ± 


.0021 


1.2073 ± 


.0016 


1.1954 ± 


.0004 


5. 


91.30 


1.1819 ± 


.0003 


1.1790 ± .0010 


1.1879 ± 


.0021 


1.1860 ± 


.0016 


1.1753 ± 


.0004 


6. 


91.95 


.8924 ± 


0003 


.8909 ± .0007 


.9058 ± 


0016 


.8965 ± 


0012 


.8899 ± 


0003 


7. 


93.00 


.4964 ± 


0003 


.4954 ± .0005 


.5004 ± 


0009 


.4976 ± 


0008 


.4971 ± 


0001 


8. 


93.70 


.3717 ± 


0003 


.3704 ± .0004 


.3690 ± 


0006 


.3720 ± 


0007 


.3730 ± 


0001 
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Figure 3: Monte Carlo results for the CALO ES, for two values (10° and 25°) of acoliinearity cut. 
Center of mass energies (in GeV) close to Z peak. In the plots, the cross section ctref from TOPAZO 
is used as a reference cross section. Cross sections in nb. 
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No. 


BHWIDE 


TOPAZO 


BHAGENE3 


UNIBAB 


SABSPV 


BHAGEN95 




(a) CALO acoLa^c = 10° 


i. 


35.257 ± .040 


35.455 ± .024 


34.690 ± .210 


34.498 ±.157 


35.740 ± .080 


35.800 ±.019 


2. 


29.899 ± .034 


30.024 ± .020 


28.780 ± .170 


29.189 ±.134 


30.270 ± .070 


30.296 ±.016 


3. 


25.593 ± .029 


25.738 ± .015 


24.690 ± .150 


24.976 ±.115 


25.960 ± .060 


25.958 ±.014 




(b) CALO acoLa,, = 25° 


1. 


39.741 ± .049 


40.487 ± .025 


39.170 ±.280 


39.521 ±.158 


40.240 ± .100 


40.463 ± .021 


2. 


33.698 ± .042 


34.336 ± .017 


32.400 ± .190 


33.512 ±.135 


34.100 ± .080 


34.287 ±.018 


3. 


28.929 ± .036 


29.460 ± .013 


27.840 ± .160 


28.710 ±.116 


29.280 ± .070 


29.409 ±.015 
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Figure 4: Monte Carlo results for the CALO ES, for two values (10° and 25°) of acollinearity cut. 
Center of mass energies close to VF-pair production threshold {Ecm- 1- 175 GeV, 2. 190 GeV, 
3. 205 GeV). In the plots, the cross section ctref from BHWIDE is used as a reference cross section. 
Cross sections in pb. 
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